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We study the influence of a temperature-dependent shear viscosity over entropy density ratio n/s, 
different shear relaxation times t„ , as well as different initial conditions on the transverse momentum 
spectra of charged hadrons and identified particles. We investigate the azimuthal flow asymmetries 
as a function of both collision energy and centrality. The elliptic flow coefficient turns out to be 
dominated by the hadronic viscosity at RH1C energies. Only at higher collision energies the impact 
of the viscosity in the QGP phase is visible in the flow asymmetries. Nevertheless, the shear viscosity 
near the QCD transition region has the largest impact on the collective flow of the system. We also 
find that the centrality dependence of the elliptic flow is sensitive to the temperature dependence 
of n/s. 
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I. INTRODUCTION 

Determining the properties of the quark-gluon plasma 
(QGP) is nowadays one of the most important goals in 
high-energy nuclear physics. For a system of weakly 
interacting particles reliable results can be obtained 
from first-principle quantum field-theoretical calcula- 
tions. Unfortunately, for strongly interacting matter 
these tools provide only a limited amount of information. 
It is, however, possible to calculate the thermodynamical 
properties of such matter numerically from the theory of 
strong interactions, quantum chromodynamics (QCD). 
These lattice QCD calculations show that if the temper- 
ature is sufficiently high, the matter undergoes a tran- 
sition from a confined phase where the relevant degrees 
of freedom are hadrons, to a deconfined phase where the 
degrees of freedom are quarks and gluons, the so-called 
QCD transition [l(. 

In recent years, experiments at the Relativistic Heavy- 
Ion Collider (RHIC) at Brookhaven National Labora- 
tory and the Large Hadron Collider (LHC) at CERN 
have provided a wealth of data from which one could in 
principle obtain information about the QGP. However, 
to compare these data with lattice QCD results is not 
straightforward. So far, lattice calculations have pro- 
vided reliable results for static thermodynamical prop- 
erties of QCD matter, e.g. the equation of state (EoS). 
The system created in heavy-ion collisions is, however, 
not static but dynamical, because it expands and cools in 
a very short time span of order 10 -23 seconds. Obviously, 
in order to be able to properly interpret the experimental 
results and infer the properties of QCD matter, we also 
need a good understanding of the dynamics of heavy-ion 
collisions. 

Fluid dynamics is one of the most commonly used 
frameworks to describe the space-time evolution of the 
created fireball, because the complicated microscopic dy- 



namics of the matter is encoded in only a few macroscopic 
parameters like the EoS and the transport coefficients. 

Currently, fluid-dynamical models give a reasonably 
good quantitative description of transverse momentum 
spectra of hadrons and their centrality dependence [3-0] . 
So far, most calculations assume that the shear viscosity 
to entropy density ratio rj/s is constant, and they show 
that, in order to describe the azimuthal asymmetries of 
the spectra, e.g. the elliptic flow coefficient V2, this con- 
stant must be very small, of order 0.1. However, for real 
physical systems, rj/s depends (at least) on the temper- 
ature [8]. A constant value of rj/s can only be justified 
as an average over the space-time evolution of the sys- 
tem. It is not clear how this average is related to the 
temperature dependence of r//s. 

In previous work [9|, [HJ, we have studied the conse- 
quences of relaxing the assumption of a constant rj/s. 
We found that the relevant temperature region where the 
shear viscosity affects the elliptic flow most varies with 
the collision energy. At RHIC the most relevant region 
is around and below the QCD transition temperature, 
while for higher collision energies the temperature region 
above the transition becomes more and more important. 
In this work we shall extend our previous study and pro- 
vide a more detailed picture of the temperature regions 
that affect elliptic flow as well as higher harmonics at a 
given collision energy. 

This paper is organized in the following way. In Sec. 
II, we describe our fluid-dynamical framework and its nu- 
merical implementation. In Sec. Ill, we specify the EoS, 
the transport coefficients, and the initialization. Sections 
IV and V contain a detailed compilation of our results, 
some of which were already shown in Refs. [3, [l(|. We 
present the transverse momentum spectra and the ellip- 
tic flow of hadrons at various centralities with different 
parameterizations of rj/s as function of temperature. We 
also study the impact of different initial conditions and 



of the choice of the relaxation time for the shear-stress 
tensor. In Sec. VI, we investigate evolution of the elliptic 
flow in more detail and, in Sec. VII, find the tempera- 
ture regions where V2 and «4 are most sensitive to the 
value of r\j s. Finally, we summarize our results and give 
some conclusions. We use natural units h = c = k = 1 
throughout the paper. 



II. FLUID DYNAMICS 

A. Formalism 

In order to describe the evolution of a system on length 
scales much larger than a typical microscopic scale, for 
instance the mean-free path, it is sufficient to characterize 
the state of matter by a few macroscopic fields, namely 
the energy-momentum tensor T^ v and, possibly, some 
charge currents N£. Fluid dynamics is equivalent to the 
local conservation laws for these fields, 



Q^V = Q 



W = 0. 



(1) 



In the absence of conserved charges and bulk viscosity, 
the energy-momentum tensor T^ v can be decomposed as 



rpfj.l> 



PA^ 



nu 



(2) 



where u^ = T MI/ w„/e is the fluid four-velocity, e is the 
energy density in the local rest frame of the fluid, i.e., in 



the frame where 



(1,0,0,0), and P is the thermo- 



dynamic pressure. The shear-stress tensor is defined as 
^iiv _ rp(^) ^ wnere t ne angular brackets () denote the 
symmetric and traceless part of the tensor orthogonal to 
the fluid velocity. With the (+,—,—,—) convention for 
the metric tensor g^ v , the projector A^ v = g^ u — u^u v . 

If the system is sufficiently close to local thermody- 
namical equilibrium, the energy-momentum conservation 
equations can be closed by providing the EoS, P{T), the 
equations determining ir 1 "' , and the transport coefficients 
entering these equations, e.g. the shear viscosity rj{T). 
The EoS P(T) and the shear viscosity r](T) can in prin- 
ciple be computed by integrating out the dynamics on 
microscopic length scales. 

While the conservation laws are exact for any system, 
the equations determining the shear-stress tensor require 
certain approximations, so that the only variables en- 
tering the equations of motion are those that appear in 
the energy- momentum tensor, namely e, u M , and it' 11 ' . 
In the so-called relativistic Navier-Stokes approximation, 
the shear-stress tensor is directly proportional to the gra- 
dients of the four- velocity, 



iV = 2r i (j>"' = 2r,d^i 



"> 



(3) 



We note that in this approximation the shear-stress ten- 
sor is not an independent dynamical variable. 

Unfortunately, this approximation results in parabolic 
equations of motion, and subsequently the signal speed is 



not limited in this theory. In relativistic fluid dynamics 
this violation of causality leads to the existence of linearly 
unstable modes, which make relativistic Navier-Stokes 
(NS) theory useless for practical applications [llllliq . 

A commonly used approach that cures these instabil- 
ity and acausality problems is Israel- Stewart (IS) theory 
[l3j ]. In this approach the shear-stress tensor, the heat 
flow and bulk viscous pressure are introduced as inde- 
pendent dynamical variables and fulfill coupled, so-called 
relaxation-type differential equations of motion. Assum- 
ing vanishing heat-flow and bulk viscosity, the relaxation 
equation for the shear-stress tensor can be written as [1 51 ] , 



TttIT 



(jiv) 



+ it^ = 277*7'"' + Xm^e + A 2 CT^ Q 7T^ Q 



(4) 



where A — u^d^A denotes the comoving derivative of 
A and 9 = d^u^ is the expansion scalar. The shear- 
relaxation time TV is the slowest time scale of the un- 
derlying microscopic theory [14| . Formally, IS theory 
can be derived by neglecting all faster microscopic time 
scales [15| . Like t„. , the coefficients Aj can in principle be 
calculated from the underlying microscopic theory, i.e., 
in our case QCD. Unfortunately, for QCD the transport 
coefficients appearing in Eq. ([2]) are still largely unknown. 
For the sake of simplicity, in this work we use Ai = —4/3, 
obtained from the Boltzinarm equation for a massless 
gas [13], and A 2 = A3 = A4 = 0. The shear- relaxation 
time and the shear viscosity are left as free parameters. 
Instead of the full (3+l)-dimensional treatment we 
consider a simplified evolution where the expansion in 
the z-direction is described by boost-invariant scaling 
flow [HI, i.e., the longitudinal velocity is given by v z = 
z/t, and the scalar densities are independent of the space- 
time rapidity r] 8 = \ log I j^ 2 - ) . Here, t is the time 
measured in laboratory coordinates. In this approxima- 
tion the full evolution depends only on the coordinates 
(r,x,y), where x and y are the transverse coordinates 
and t — \/t 2 — z 2 is the longitudinal proper time. 



B. Numerical implementation 

Once the initial values of the components of the energy- 
momentum tensor are specified at a given initial time To, 
the space-time evolution of the system is obtained by 
solving the conservation laws ([T]) together with the IS 
equations dU). 

The conservation laws are solved using the algorithm 
developed in Refs. 11711 and generalized to more than one 
dimension in Ref. [ljj. This method, known as SHASTA 
for "SHarp and Smooth Transport Algorithm", solves 
equations of the type 



dtU + di(viU) = S(t,x) 



(5) 



where U — U(t,x) is for example T 00 , T *, . . ., Vi is the 
ith component of three-velocit y, a nd S(t, x) is a source 
term, for more details see Ref. [19j . 



We can further stabilize SHASTA by letting the an- 
tidiffusion coefficient A a d which controls the amount of 
numerical diffusion to be proportional to 



1 



(k/ey + 1 



(6) 



where e is the energy density in the local rest frame, and 
k is some constant of order 10 -5 GeV/fm 3 . In this way, 
Aad goes smoothly to zero near the boundaries of the 
grid, i.e., we increase the amount of numerical diffusion in 
that region. We have checked that this neither affects the 
solution nor produces more entropy inside the decoupling 
surface. 

The relaxation equation (j4j could also be solved using 
SHASTA. However, we noticed that solving it by replac- 
ing the spatial gradient at grid point i on the left-hand 
side of Eq. (|4]) by a centered second-order difference, 



d x Ui 



U iA 



Vi. 



2A.t 



(7) 



where U = tt^ , yields a more stable algorithm. Time 
derivatives in the source terms are simply taken as first- 
order backward differences. Like in SHASTA, all spatial 
gradients in the source terms are discretizised according 
to Eq. 0. 



III. PARAMETERS 

A. Equation of State 

As EoS we use the recent s95p-PCE-vf parameteriza- 
tion of lattice QCD results [23J. In this parameteriza- 
tion, the high-temperature part is matched to recent re- 
sults of the hotQCD collaboration [2J, |25| and smoothly 
connected to the low-temperature part described as a 
hadron resonance gas. All hadrons listed in Ref. [26| up 
to a mass of 2 GeV are included in the hadronic part of 
the EoS. The system is assumed to chemically freeze-out 
at T c hcm = 150 MeV. Below this temperature the EoS 
is constructed according to Refs. |27H29{ . This construc- 
tion assumes that the evolution below T c h e m is isentropic. 
Strictly speaking this is not the case in viscous hydrody- 
namics since dissipation causes an increase in entropy. 
However, we have checked that in our calculations the 
viscous entropy production from all fluid cells with tem- 
peratures below T c hcm = 150 MeV is less than 1% of the 
initial entropy, whereas the entropy production during 
the entire evolution ranges from 3 — 14 %, depending on 
the collision energy and the 77/s parameterization. 



B. Transport coefficients 



C. Freeze-out 

We assume that freeze-out, i.e., the transition from the 
fluid-dynamical system to free-streaming particles hap- 
pens on a hypersurface of constant temperature. Unless 
otherwise stated, we assume that the freeze-out temper- 
ature is Tdcc = 100 MeV. We include all 2- and 3-particle 
decays of hadronic resonances according to Ref. [2fJ . 

The transverse momentum distribution of hadrons 
is calculated using the Cooper- Frye description |21j . 
For the final spectra we need to know the lo- 
cal single-particle momentum distribution functions of 
hadrons on the freeze-out surface. Here, we em- 
ploy the widely used I4-moment ansatz where the 
correction to the local-equilibrium distribution /qj = 
{exp [{u^pl — fii) jT\ ±1} of a hadron of species i with 
four-momentum p% is given by [22| 



Sfi = /, 



0/ 



T 2 (e + P) ■ 



(8) 



We note that this functional form for Sf is merely an 
ansatz. If dissipative fluid dynamics is derived from the 
Boltzmann equation without assuming the 14-moment 
approximation, the full expansion of Sf contains an infi- 
nite number of terms, for details see Ref. [15|. The effect 
of this will be studied in a future work. 



The temperature-dependent shear viscosity is 
parametrized as follows. In all cases, we take the 
minimum of 77/s to be at T tr = 180 MeV. Unless other- 
wise stated, the value of r\j ' s at the minimum is assumed 
to equal the lower bound 77/s = 0.08 conjectured in the 
framework of the AdS/CFT correspondence [311 ] . 
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FIG. 1. (Color online) Different parameterizations of rj/s as a 
function of temperature. The (LH-LQ) line is shifted down- 
wards and the (HH-HQ) line upwards for better visibility. 

The parameterization of the hadronic viscosity is based 
on Ref. [32| where the authors consider a hadron reso- 
nance gas with additional Hagedorn states. In practice, 
we use a temperature dependence of 77/s of the following 



functional form [9|, (33| , 



= 0.681-0.0594— - 0.544 f— I . (<)] 
HRG Ttt \r tr ' 



At T = 100 McV this coincides with the rj/s value given 
in Ref. [32j . and decreases smoothly to the minimum 
value rj/s — 0.08 at T tr . We note that many authors ob- 
tain considerably larger values for the shear viscosity of 
hadronic matter, see e.g. Refs. [3J|. Our motivation here 
is to illustrate the effects of hadronic viscosity rather than 
to use a parameterization that is as realistic as possible. 
We shall see that even this low rj/s leads to considerable 
effects for hadronic observables in Au + Au collisions at 
RHIC. We further note that, since we are considering 
a chemically frozen hadron resonance gas below T c h Gm , 
while in Ref. |32J chemical equilibrium is assumed at all 
temperatures, the entropy densities, and therefore the 
values of rj, differ between the two calculations at a given 
value of T < T c hem- 

The high-temperature QGP viscosity is parametrized 
according to lattice QCD results [35| in such a way that 
it connects to the minimum of rj/s at T tr . The functional 
form used is 



QGP 



= -0.289 



T 
0.288 — 



0.0818 



h 



We take the following four parameterizations of the 
shear viscosity: 

• (LH-LQ) rj/s = 0.08 for all temperatures, 

• (LH-HQ) rj/s = 0.08 in the hadron gas, and above 
T = 180 MeV rj/s increases according to Eq. (TTU|) , 

• (HH-LQ) below T = 180 MeV, rj/s is given by Eq. 
©, and above we set rj/s = 0.08, 

• (HH-HQ) we use Eqs. © and (JTUJ) for the hadron 
gas and the QGP, respectively. 

These parameterizations are shown in Fig. [1] Besides 
these four cases we also study the effect of varying the 
value of the minimum of rj/s, see Sees. IVl and IVIII 

In order to complete the description, we also need to 
specify the relaxation time. In this work we use a func- 
tional form suggested by kinetic theory, 



(11) 



P 



where c T is a constant. Causality requires that c T > 
2 [I3|. Unless otherwise stated, we shall use the value 
c T = 5 which coincides with the value obtained from the 
Boltzmann equation in the 14-momcnt approximation for 
a massless gas of classical particles [3g|. The relaxation 
times corresponding to the parameterizations above are 
shown in Fig. [2] The effect of varying the relaxation time 
separately from r/ is also studied in Sec. [Vj 
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FIG. 2. (Color online) Relaxation times corresponding to the 
different parameterizations of rj/s, for c T = 5. The (LH-LQ) 
line is shifted downwards and the (HH-HQ) line upwards for 
better visibility. 

C. Initial state 

We still need to specify the initial state at some proper 
time To- For a boost-invariant system it is sufficient to 
provide the components of the energy-momentum tensor 
in the transverse plane at z = 0, i.e., tj s = 0. Within 
our approximations these are the local energy density, 
the initial transverse velocity, and the three independent 
components of the shear-stress tensor. Here, we will as- 
sume that the initial transverse velocity is zero and, un- 
less otherwise stated, the initial shear-stress tensor is also 
assumed to be zero. 

For the initial time we choose To = 1 fm. The energy 
density e(ro,x,y) is based on the optical Glauber model 
by assuming that the energy density is a function of the 
density of binary nucleon-nucleon collisions tt-bc? o r the 
density of wounded nucleons n-yvN, or both, 



e(r , x, y) = C e /(n B c, ™wn). 



(12) 



The overall normalization, C e , is fixed in order to re- 
produce the observed multiplicities in the most central 
^Jsnn = 200 GeV Au+Au collisions at RHIC, and in 
V^iv = 2.76 GeV Pb+Pb collisions at LHC. 

The centrality dependence of the multiplicity is repro- 
duced in this work in two different ways: 

• BCfit: choosing / to be a polynomial in Ubc, 

f(mc) =^BC +Ci7l| c + C 2 n| c . (13) 

• GLmix: using a superposition of nee and n-yvN, 

/("bc^wn) = din BC + (1 -di)nwN- (14) 

Here, the coefficient ci is introduced in order to guarantee 
that the parameterizations are monotonically increasing 
with increasing binary-collision or wounded-nucleon den- 
sity. This ensures that the highest energy density is in 
the center of the system, i.e., at x = y = 0. 
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c 2 [fm- 4 ] 


di 


T max [MeV] 


200 


-0.032 


0.00035 


0.1 


313 


2760 


-0.01 


0.0001 


0.7 


430 


5500 








1.0 


504 



TABLE I. Initialization parameters for different collision en- 
ergies. The maximum temperature T ma x is given for the BCfit 
initialization with the (LH-LQ) parameterization of r]/s. For 
the other initializations T max differs less than 5%. 

For a given impact parameter, the optical Glauber 
model yields a different number of participants and differ- 
ent centrality classes than the Monte Carlo Glauber mod- 
els commonly used by the experimental collaborations. 
Using the optical Glauber model, we can either choose 
to reproduce the multiplicity as a function of the num- 
ber of participants or as a function of centrality classes. 
In general, this leads to different coefficients Ci and d\. 
Here, we choose to determine the initial conditions by 
requiring that the centrality dependence of the charged 
particle multiplicity as a function of the number of partic- 
ipants J37|,|38| is reproduced. We have checked that, if we 
determine the centrality dependence by matching to the 
centrality classes given by the optical Glauber model, the 
elliptic flow is more suppressed in central and enhanced 
in peripheral collisions at RHIC energies, while at LHC 
energies it remains practically unchanged. In order to be 
fully consistent with the experimental determination of 
the centrality classes, one would need to generate fluctu- 
ating initial conditions via a Monte Carlo Glauber model, 
see e.g. Refs. HHU- 

For ^/smn = 5.5 TeV Pb+Pb collisions we use the mul- 
tiplicity in the most central collisions as predicted by the 
EKRT model [4J| . In this case the centrality dependence 
is assumed to follow binary scaling, i.e., c\ = c 2 = 
in Eq. (|13p . All initialization parameters are shown in 
Table H 

Different parameterizations of r\j s lead to different en- 
tropy production and therefore different final multiplic- 
ity, even if the initial state is kept the same. This is es- 
pecially true for different parameterizations of the high- 
temperature shear viscosity, since most of the ent ropy 
is produced during the early stages of the collision (42j | . 
We compensate this using different overall normalizations 
e.g. between the (HH-LQ) and (HH-HQ) parameteriza- 
tions. Entropy production during the hadronic evolution 
is small and not compensated. The centrality dependence 
of the entropy production is also different for different r)/s 
parameterizations. Since it leads to at most a 5% differ- 
ence in the final multiplicities and is hardly visible in the 
results, it is not corrected here. 



IV. RESULTS AND COMPARISON WITH 
EXPERIMENTAL DATA 

In this section we use the initializations and parame- 
terizations of r//s given above, and compare the results 



with experimental data from RHIC and LHC. 



Transverse momentum spectra and elliptic flow 
at RHIC 



In Fig. [3] we show the pr-spectra of pions for different 
centrality classes for RHIC ^/snn = 200 GeV Au+Au 
collisions and compare them with PHENIX data [37[ . We 
only show results using the BCfit initialization; those for 
the GLmix initialization are very similar. The freeze-out 
temperature is chosen as T^ec = 100 MeV. This choice 
reproduces the slopes of the pT-spectra quite well. 

Once we correct the normalization of the initial energy 
density profile for different entropy production, the slopes 
of the pT-spectra are practically unaffected by the n/s pa- 
rameterizations. We note that in our earlier work [9| this 
correction was not made, and the different rj/s param- 
eterizations lead not only to different multiplicities but 
also to different slopes for the pT-spectra. This effect was 
even more pronounced at LHC than at RHIC, due to an 
increase in entropy production caused by larger gradients 
appearing with an earlier initialization time tq — 0.6 fm. 
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FIG. 3. (Color online) Pion spectra at RHIC, with BCfit 
initialization. 

The kaon spectra are shown in Fig. |4] and the proton 
spectra in Fig. [5] with the BCfit initialization. Both are 
compared with PHENIX data 37[. Because we do not 
consider net-baryon number in our calculations, the pro- 
ton and anti-proton spectra are identical. For this rea- 
son we show both the proton and the anti-proton data in 
Fig. El 

For both kaons and protons the calculated spectra are 
slightly more curved than the data and they also lie above 
the data. As for the pions, the slopes of the spectra are 



practically independent of the 77/ s parameterization. 




FIG. 4. (Color online) Kaon spectra at RHIC, with BCfit 
initialization. 



at RHIC. This holds for all centrality classes. The sup- 
pression of the elliptic flow due to the hadronic viscos- 
ity is even more enhanced in more peripheral collisions. 
Note that with the BCfit initialization, the elliptic flow 
in the most central collision class is reproduced by the 
parameterizations with a large hadronic viscosity, while 
with the GLmix initialization the elliptic flow in the same 
centrality class is better described by taking a constant 
rj/s = 0.08. However, with the latter choice the elliptic 
flow tends to be overestimated in more peripheral col- 
lisions. On the other hand, the temperature-dependent 
hadronic rj/s gives the centrality dependence correctly 
down to the 30 — 40 % centrality class. In even more 
peripheral collisions a large hadronic viscosity tends to 
suppress the elliptic flow too much. 

Figure [8] shows V2(pt) for protons with the BCfit ini- 
tialization compared to the two-particle cumulant data 
from the STAR collaboration J44j]. The protons show 
qualitatively the same response to the different rj/s pa- 
rameterizations as all charged hadrons, i.e., V2(pt) de- 
pends strongly on the hadronic viscosity, but is almost 
independent of the high-temperature rj/s. Since we use a 
smooth initialization, with no initial-state fluctuations 
included, quantitative comparisons with two- or four- 
particle cumulant data are not straightforward. 
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FIG. 5. (Color online) Proton spectra at RHIC, with BCfit 
initialization. 

Figure[6]shows the pr-differential elliptic flow V2(pt) of 
charged hadrons for different centrality classes using the 
BCfit initialization. Similarly, Fig. [7] shows the elliptic 
flow for the GLmix initialization. The calculations are 
compared with the four-particle cumulant data from the 
STAR collaboration 0. 

As was observed in Ref. [9| , the differential elliptic flow 
is largely independent of the high-temperature rj/s pa- 
rameterization, but highly sensitive on the hadronic rj/s 
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FIG. 6. (Color online) Charged hadron «2(pr) at RHIC, with 
BCfit initialization. 



B. Transverse momentum spectra and elliptic flow 
at LHC 

Transverse momentum spectra of charged hadrons in 
most central Pb+Pb collisions with ^/snn = 2.76 TeV 
at LHC are shown in Fig. [S] At LHC, both initializa- 
tions BCfit and GLmix give very similar results for both 
elliptic flow and the spectra, because the contribution 




1 2 

Pt [QeV] 



FIG. 7. (Color online) Charged hadron «2(pt) at RHIC, with 
GLmix initialization. 



£ 0.10 



0.05 



0.15 



it 0.10 




RHIC 200 AGeV 
H >75 H 



30-40 % 





2 



Pt [GeV] 



Pt [GeV] 



FIG. 8. (Color online) Proton u 2 (pt) at RHIC, with 5C/ii 
initialization. 



from binary collisions is large, of order ~ 70 %, see Ta- 
bic HI Therefore, we show only results with the BCfit ini- 
tialization; these are compared to data from the ALICE 
collaboration [45J. The calculated spectra are somewhat 
flatter than the data. Here, we have used the same decou- 
pling temperature as at RHIC, i.e., Td cc = 100 MeV. We 
could improve the agreement with the data by decoupling 
at even lower temperature than at RHIC. Another way 
to improve the agreement is choosing a larger chemical 
freeze-out temperature. This would give steeper spectra, 
but the proton multiplicity at RHIC would then be over- 
estimated. However, we have tested that the dependence 
of the spectra and the elliptic flow on rj/s is unchanged 



by these details. 

As was the case at RHIC, at LHC the slopes of the 
spectra are practically independent of the rj/s parame- 
terization. We note that here we have used the initial- 
ization time to = 1.0 fm, i.e., the same as at RHIC. In 
Ref. [9] we observed a quite visible correlation between 
the shear viscosity and the spectral slopes. Here, the later 
initialization time and the fact that we now compensate 
for the entropy production between different rj/s param- 
eterizations almost completely removes this correlation. 
However, the earlier the evolution starts, the more the 
viscosity will affect the slopes. 
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FIG. 9. (Color online) Charged hadron spectra at LHC, with 
BCfit initialization. 

The pr-differential elliptic flow for all charged hadrons 
is shown in Fig. QJj] and for protons in Fig. [TTJ The 
charged hadron elliptic flow is compared with ALICE 
four-particle cumulant data |46|. We can see that in the 
10 — 20 % centrality class, changing the hadronic rj/s or 
changing the high-temperature rj/s has quite a similar 
impact on the elliptic flow, e.g. the difference between 
the LH-LQ and LH-HQ and between the LH-LQ and 
HH-LQ curves is nearly the same. However, the more pe- 
ripheral the collision is, the more the viscous suppression 
is dominated by the hadronic rj/s. This is confirmed by 
comparing the grouping of the flow curves in the 40 — 50 
% centrality class at LHC with that at RHIC, cf. Figs. [6] 
andflOl As was the case in Au+Au collisions at RHIC, 
also here the grouping of the curves for the protons is 
similar to that of all charged hadrons, cf. Fig. [TTJ 

Note that, within our set-up, the best agreement with 
the ALICE data is obtained with the HH-HQ parameter- 
ization, i.e., with a temperature-dependent rj/s in both 
hadronic and high-temperature phases. However, in the 
low-pT region our calculations systematically underesti- 
mate the elliptic flow in all centrality classes. As was the 
case with the pr-spectrum, decoupling at a lower tern- 



perature and choosing a higher chemical freeze-out tem- 
perature would improve the agreement, without changing 
the grouping of the elliptic flow curves with the rj/s pa- 
rameterizations. 
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FIG. 12. (Color online) Charged hadron V2(pt) 
A TeV, with BC initialization. 



at LHC 5.5 



FIG. 10. (Color online) Charged hadron V2(j>t) at LHC, with 
BCfit initialization. 




FIG. 11. (Color online) Proton v 2 {p T ) at LHC, with BCfit 
initialization. 

In Fig. [T2] we show the p^-differential elliptic flow for 
y^sJfN = 5.5 TeV Pb+Pb collisions. In this case the 
viscous suppression of V2(pt) is dominated by the high- 
temperature rj/s in central collisions, while peripheral 
collisions resemble more the lower-energy central colli- 
sions at LHC, i.e., both hadronic and high-temperature 
viscosity contribute similarly to the suppression. Fur- 
thermore, the higher the px, the more the hadronic 
viscosity contributes to the suppression. This happens 
mainly because 6f increases with both viscosity and pr- 



V. EFFECTS OF SHEAR INITIALIZATION, 
MINIMUM OF n/a AND RELAXATION TIME 

One of the main results of Ref. [9j is that, at RHIC, 
the high-temperature shear viscosity has very little effect 
on the elliptic flow. In this section we elaborate more 
on this analysis, and explicitly show that this statement 
holds for an out-of-equilibrium initialization of the shear- 
stress tensor as well. We also study the effect of varying 
the relaxation time. 
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FIG. 13. (Color online) Charged hadron v 2 (p T ) at RHIC, 
with BCfit and NS initialization. 

Figure [T3] shows the elliptic flow of charged hadrons in 
the 20 — 30 % centrality class at RHIC. Instead of set- 
ting tt^ to zero initially, here the so-called Navier-Stokes 



(NS) initialization where the initial values of the shear- 
stress tensor are given by the first-order, asymptotic solu- 
tion of IS theory, Eq. ([3]). For all rj/s parameterizations, 
the NS initialization increases the entropy production (up 
to 30 %), especially for the parameterizations with a large 
high-temperature viscosity. This is corrected by adjust- 
ing the initial energy density to produce approximately 
the same final multiplicity. Although for the parame- 
terizations with a large hadronic rj/s the different shear 
initializations give slightly different v 2 (pt) curves, the 
grouping of these curves remains intact. We emphasize 
that the NS initialization gives very different initial con- 
ditions for each viscosity parameterization. If we use the 
same non-zero initial shear stress, e.g. ir^ 1 " = const, x a^ v , 
for each parameterization, the resulting v 2 {j>t) curves in 
each group in Fig. Q2] would be even closer to each other. 
The NS initialization with a constant rj/s — 0.08 has 
a relatively short relaxation time, see Fig. [2] Hence for 
T n <C To the NS initialization is not a completely unrealis- 
tic assumption for the initial values of it 111 ' . However, for 
larger values of rj/s the relaxation times are considerably 
larger, r n > To, and there is no reason to assume that 
the asymptotic solution could have been reached already 
at very early times. 
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FIG. 14. (Color online) Parameterizations of rj/s as a function 
of temperature. The (HH-HQ) line the same as in Fig. [TJ 

So far we have changed the shear- viscosity parameteri- 
zation by keeping the minimum fixed. In Fig. [14] we show 
the original HH-HQ parameterization and one where rj/s 
around the minimum is twice as large. Figure [15] shows 
three v 2 (pt) curves for Au+Au collisions at RHIC: one 
with the original HH-HQ parameterization, one with the 
larger minimum value of rj/s, and the last one with the 
same large minimum value of rj/s, but with a larger re- 
laxation time, i.e., the constant in the relaxation time 
formula (TTTj) is c T — 10 instead of c T = 5. We note that 
even a relatively small change in the rj/s parameteriza- 
tion near the minimum produces quite a visible change 
in v 2 (pt). At RHIC, this change can be almost com- 
pletely compensated by adjusting the relaxation time. 
This shows that in small, rapidly expanding systems like 
the one formed in heavy-ion collisions, transient effects 



have considerable influence on the evolution. In other 
words, the relaxation time cannot be merely considered 
as a way to regularize the unstable Navier-Stokes theory: 
it has real physical effects that cannot be completely dis- 
tinguished from the effects of rj/s. In ^/smn = 2.76 TeV 
Pb+Pb collisions at LHC, the effect of changing the min- 
imum or the relaxation time is practically the same. 
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FIG. 15. (Color online) Charged hadron V2(pr) at RHIC, 
with BCfit initialization and for different minima of rj/s and 
relaxation times. 



VI. TIME EVOLUTION OF THE ELLIPTIC 
FLOW 

One way to probe the effects of shear viscosity on the 
elliptic flow is to calculate the time evolution of the lat- 
ter. Typically this is done by calculating the so-called 
momentum-space anisotropy from the energy-momentum 
tensor, 



irpxx _ j^yy\ 
jrpxx _|_ rpyy) ' 



(15) 



where the (• • ■ ) denotes the average over the transverse 
plane. The problem is, however, that one cannot make a 
direct connection of s p to the actual value of v 2 obtained 
from the decoupling procedure. Also, this way of study- 
ing the time evolution does not take into account that, 
at fixed time, part of the matter is already decoupled, 
i.e., the average over the transverse plane includes also 
matter that is outside the decoupling surface. 

To overcome these two shortcomings of e p , we in- 
stead calculate the v 2 of pions from a constant-time 
hypersurface that is connected smoothly to a constant- 
temperature hypersurface at the edge of the fireball, see 
Fig.HHfor examples of such hypersurfaces. Although, the 
pions do not exist as real particles before hadronization, 
the advantage is that the final v 2 we obtain matches the 
one of thermal pions from the full decoupling calculation. 

Figure 1171 shows the time evolution of v 2 in Au+Au 
collisions at RHIC, in ^/s~Wn =2.76 TeV Pb+Pb colli- 
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FIG. 16. (Color online) Constant-temperature hypersurfaces 
at decoupling (Td cc = 100 MeV) , chemical freeze-out (T c hem = 
150 MeV), and at the minimum of r)/s (T tr = 180 MeV) at 
different collision energies. Also, examples of surfaces that are 
used in the calculation of the time evolution of V2 are shown 
(dotted lines). 



sions at LHC, and in ^/snn — 5.5 TeV Pb+Pb collisions 
at LHC. In all cases, the evolution is calculated in the 
20 — 30 % centrality class. These results confirm our 
earlier conjecture: at RHIC, the different rj/s parameter- 
izations create very little difference in the elliptic flow in 
the early stages of the collision, while at later stages the 
suppression due to the hadronic viscosity takes over and 
groups the V2 curves according to the hadronic viscosity. 
At the intermediate LHC energy the impact of the QGP 
viscosity is larger, and the final V2 still has a memory 



FIG. 17. (Color online) Time evolution of V2 at different 
collision energies. 



of this difference. The hadronic viscosity has a similar 
impact on V2 as the QGP viscosity. At the highest LHC 
energy the hadronic suppression is small and the effect 
of the QGP viscosity clearly dominates the grouping of 
the i>2 curves. Interestingly, both LHC evolutions show 
an increase of v-i around r = 10 fm/c. This is when the 
system is going through the chemical decoupling stage. 
In the chemically frozen system V2 tends to increase more 
rapidly than in chemical equilibrium [23. |30I| - At RHIC, 
the chemical decoupling happens earlier, and also the 
hadronic suppression is stronger, and the increase in V2 
is washed out. 
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VII. PROBING THE EFFECTS OF A 
TEMPERATURE-DEPENDENT rj/s ON THE t>„'S 

In this section, we try to probe the effects of a 
temperature-dependent 77/s on the azimuthal asymme- 
tries in a more detailed way. To this end, we introduce 
a modified 77/s. Our baseline is a constant rj/s\ c = 0.08 
that we then modify near some temperature Ti according 
to 



V = V 

s s 



1 + 2 exp 



\T-TA-ST 



1 



(16) 

where the parameters are taken to be ST = 10 MeV and 
A = 1.5 MeV. One example of this 77/s parameteriza- 
tion is shown in Fig. [18] We note that, although we 
use smooth initial conditions from the optical Glauber 
model, we still get non-zero v n for all even n. Although 
these are much smaller than the ones obtained with the 
fluctuations included, we can still probe the effects of vis- 
cosity on these coefficients. By changing the temperature 
Ti and comparing the simulations with a constant 77/s we 
can find the temperature regions where 1*2 or 7J4 are most 
sensitive to changes of 77/s at different collision energies. 
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FIG. 18. (Color online) Shear viscosity with a modified tem- 
perature dependence. 

Figure [TO] shows the results for vi and 7J4 in the 
20 - 30 % centrality class for RHIC and for both LHC 
energies considered earlier. We plot the relative differ- 
ence Sv n /v n , where Sv n = v n (r)/s(T)) - tj„(t7/s| c ). Each 
point in the figure corresponds to a different calculation, 
with a different value of T in Eq. (TIB"]) . Similarly, Fig. l20l 
shows the same result, but without the Sf contribution 
to the freeze-out. 

The viscosity can affect v„ in two ways: by chang- 
ing the space-time evolution of the integrated quantities 
like the energy density, or by changing the local particle- 
distribution function at freeze-out. With our small base- 
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FIG. 19. (Color online) Effects of modified rj/s on V2 and V4. 



line viscosity the effect on the local distribution func- 
tion is quickly washed out during the evolution below 
the temperature Tj. Therefore, in these simulations, in 
most of the temperature points, the change in 77/s af- 
fects v n through the space-time evolution, except at the 
lowest-temperature point Ti = 110 MeV, where the peak 
in 77/s is close to the freeze-out temperature T^ ec = 100 
MeV. If we exclude the lowest temperature point in V4 at 
RHIC, we can read off from the figures that the temper- 
ature region where viscosity affects both V2 and V4 most 
is around the transition region T ~ 150 . . . 200 MeV. For 
V2 this temperature region shifts slightly towards higher 
temperatures with increasing collision energy, while for V4 
the temperature where the effect is maximal is practically 
unchanged. Other than this, the overall behavior of v 2 
and H4 is quite similar. At high temperatures, the effect 
of 77/ s increases with increasing collision energy, while at 
low temperatures the viscous suppression decreases with 
increasing collision energy, which is most notable for the 
Ti = 110 MeV point where the viscosity effects on the 
freeze-out distribution are strongest. 

For V2 we observed earlier that the suppression due to 
the hadronic viscosity practically vanishes at the highest- 
energy LHC collisions. This is again confirmed in Fig. 1191 
This is, however, not true for higher harmonics. For tj 4 
there is still a significant contribution from hadronic vis- 
cosity at the full LHC energy. In this sense, higher har- 
monics do not give direct access to the high-temperature 
viscosity, but can rather help in constraining the hadronic 
dynamics and viscosity as well as the correct form of Sf. 
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FIG. 20. (Color online) Same as Fig. QU but without the Sf 
contribution. 



This is also important since the hadronic evolution al- 
ways tends to shadow the effects of the properties of the 
high-temperature matter. 



VIII. CONCLUSIONS 

We have studied the effects of a temperature- 
dependent r\f ' s on the azimuthal asymmetries of hadron 
transverse momentum spectra. We found earlier [9j that 
the viscous suppression of the elliptic flow is dominated 
by the hadronic viscosity in ^/snn = 200 GeV Au+Au 
collisions at RHIC, while in Pb+Pb collisions at the full 
LHC energy ^/snn = 5.5 TeV the suppression is mostly 
due to the high-temperature shear viscosity. In this work 
we have supplemented these earlier studies with more de- 
tails. 

First, we found that the suppression of the elliptic flow 
due to the shear viscosity becomes more important in 
more peripheral collisions. At least in our set-up, for 
RHIC energies a temperature-dependent shear viscosity 



improves the centrality dependence of the elliptic flow 
compared to the data, similarly to what was found in the 
hybrid approach of Ref. p. With a constant rj/s = 0.08 
and with the GLmix initialization, the measured i>2(j»t) 
is reproduced in the most central collisions, but the cal- 
culations give a too large elliptic flow for peripheral colli- 
sions. However, with the BCfit initialization the elliptic 
flow in the most central collisions is reproduced with a 
temperature-dependent viscosity, and also the centrality 
dependence is reproduced down to the 30 — 40 % cen- 
trality class. Similarly, in Pb+Pb collisions at LHC both 
a temperature-dependent hadronic rj/s as well as an in- 
creasing rj/s in the high-temperature phase help in re- 
producing the centrality dependence. Although there are 
lots of uncertainties associated with the decoupling and 
the initial state, at RHIC the centrality dependence of 
V2(pt) may give access to the temperature dependence 
of rj/s in hadronic matter. 

Furthermore, we have studied the effects of a 
temperature-dependent rj/s in a more detailed way. We 
found that for a given collision energy both vi and V4 
are most sensitive to the shear viscosity near the transi- 
tion temperature, i.e., T ~ 150 — 200 MeV. For V2, this 
region moves slightly to higher temperature and widens 
with increasing collision energy, while for W4 it remains 
practically unchanged. Other than that, the dependence 
of V2 and 1)4 on rj/s is similar with increasing collision 
energy: the effect of the hadronic viscosity decreases and 
the effect of the high-temperature viscosity increases. 

For V2 the effect of Sf almost vanishes at the highest 
collision energies, but for V4 it always remains signifi- 
cant. At RHIC the Sf corrections clearly dominate U4, 
and even at the highest collision energies this effect is 
comparable to the effects due to the modified space-time 
evolution. In this sense, higher harmonics give access to 
the Sf corrections and the hadronic viscosity rather than 
the high-temperature viscosity. 
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